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ABSTRACT 

O , 

O , The convective collapse of thin magnetic flux tubes in the photospheres of 

^ . sun-like stars is investigated using realistic models of the superadiabatic upper 

. convection zone layers of these stars. The strengths of convectively stable flux 



> 



tubes are computed as a function of surface gravity and effective temperature. We 
find that while stars with Te// > 5500 K and log g> 4.0 show flux tubes highly 
evacuated of gas, and hence strong field strengths, due to convective collapse. 



^ , cooler stars exhibit flux tubes with lower field strengths. Observations reveal the 

Q , existence of field strengths close to thermal equipartition limits even in cooler 



stars, implying highly evacuated tubes, for which we suggest possible reasons. 

Subject headings: stars: atmospheres — stars: magnetic fields — instabilities- 
MHD 
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d • 1. Introduction 

^> : 

^ I An important physical quantity basic to the understanding of stellar magnetism is the 

^ [ field strength in small-scale flux tubes and spots on the photospheres. Measurements and 
interpretation of magnetic field strengths (Saar 1988, 1996a; Guenther 1997) in cool main- 
sequence stars other than the Sun have been a subject of much debate recently (Basri et al. 
1990; Safier 1999), because of the difficulties associated with the modeling of the atmospheric 
structural changes that the inhomogeneous fibril state of the magnetic field introduces (Basri 
et al. 1990). However, recent improvements in observing techniques (Saar 1996a; Donati et 
al. 1997; Johns-Krull et al. 1999) have increased the reliability of observational determination 
of stellar magnetic field strengths. The observed Zeeman broadening on cool stars is believed 
to be produced by small flux tubes that appear bright, similar to the solar magnetic network 
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and facular bright points, rather than spots, which are much fainter than quiet photospheres 
and hence contribute httle to the stellar profiles (Basri et al. 1990; Safier 1999). 

Based on measurements of Zeeman effect in cool main sequence stars Saar (1994, 1996a) 
has inferred that the magnetic field strengths are close to the thermal equipartition field 
strengths Beq—y/8Trpl at the observed levels in the atmosphere, where Pe is the gas pressure 
in the unmagnetised atmosphere, with a conclusion that the surface distribution is in the 
form of highly evacuated small flux tubes in pressure balance with the ambient atmosphere. 
Various other measurements (Saar & Linsky 1985; Basri & Marcy 1994; Johns-Krull et al. 
1999) also establish such field strengths leading to the general acceptance (Linsky 1999) that 
stellar surface field strengths scale as the square root of the surface gas pressure. We note 
that Beg, at any geometrical level in a stellar atmosphere, is the maximum possible field 
strength for flux tubes confined by gas pressure and can be easily determined (Bunte & Saar 
1993; Safier 1999) as a function of log g and Te// using atmospheric models (e.g., Kurucz 
(1993)), where g and Tg// denote the surface gravity and effective temperature. 

The widely accepted mechanism to account for the kG range field strengths in solar 
magnetic fiux tubes is the superadiabatic effect or convective collapse (CC, hereafter) (Parker 
1978; Spruit & Zweibel 1979). In this Letter, we report results from a study of this mechanism 
in the stellar context and test its efficiency as a function of log g and Te//. 

Concentration of magnetic fields into discrete fiux tubes or sheaths has long been recog- 
nised as a general consequence of the interaction between magnetic fields and cellular con- 
vection that expels the fiux into downdrafts (Parker 1963; Galloway & Weiss 1981; Proctor & 
Weiss 1982). High-resolution observations of solar surface magnetic fields confirm the opera- 
tion of such flux expulsion process. The CC, which is a consequence of thermal insulation of 
the expelled magnetic flux against convective motions and the superadiabatic thermal strat- 
iflcation of the ambient atmosphere, further intensifles the fleld to observed values (Parker 
1978). It is a convective instability, modifled by the magnetic fleld, and drives a down flow 
along the field lines of a fiux concentration leading to the evacuation of gas inside; lateral 
pressure balance ensures that a highly compressed intense field fiux tube or sheath is formed. 
This instabihty develops very rapidly, typically on a free-fall time scale, and is faster than 
other MHD instabilities present. Hence, the criterion for convective stability can be regarded 
as a test to check for the occurrence of flux tubes in stars. On the Sun, this instability has 
a typical growth time of about 2-5 minutes (Hasan 1984; Rajaguru & Hasan 2000). Other 
important instabilities are the interchange instability (Parker 1975; Piddington 1975) con- 
nected with the field topology, which has been studied in the stellar context by Bunte & Saar 
(1993), and the Kelvin- Helmholtz instability (Schussler 1979; Tsinganos 1980) connected 
with the relative motion between the field confined and the surrounding gas. 
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We use model atmospheres constructed using the ATLAS9 model atmosphere code 
(Kurucz 1993), which can extend the Kurucz (1993) atmospheres to cover deeper regions of 
the convection zone, to study the CC of thin flux tubes embedded in them. The critical field 
strengths for the tubes to be stable, which indicate the amount of evacuation that stellar 
flux tubes undergo as a result of the collapse process, are found. 

2. Convective Stability of Stellar Thin Flux Tubes 

A simple and very effective modeling of the dynamics of thin fiux tubes, including 
the important effects of stratification and compressibility, is provided by the thin fiux tube 
approximation (Roberts & Webb 1978) to the MHD equations. The thin fiux tube equa- 
tions have proved very useful in the study of small-scale solar magnetic elements and have 
been successful in explaining a variety of observed phenomena (Solanki 1993 and references 
therein). We solve the same second order differential equation derived by Spruit & Zweibel 
(1979) based on linear stability analysis of thin tube equations in the adiabatic limit, which is 
valid for arbitrary background stratifications, employing a grid of stellar model atmospheres. 
We do not reproduce the equation here, and we refer the reader to Spruit & Zweibel (1979) 
for details. The stability of a fiux tube is controlled by f3, which is the ratio of gas to mag- 
netic pressure inside the tube. Weak field (i.e, high f3) flux concentrations that are formed 
at the downdrafts of convective cells are convectively unstable and collapse until their field 
strengths reach the critical value at which the instability is quenched. Thus the instability is 
self- limiting and leads to a lower energy stable equilibrium state for flux tubes (Spruit 1979). 

2.1. Stelleir Model Atmospheres 

Models of stellar convective envelopes that extend downward to depths required to study 
the CC have been specially constructed, extending the Kurucz (1993) models to cover deeper 
layers using the ATLAS9 model atmosphere code. The convective flux is calculated using 
mixing length theory. We computed models, with solar metallicity, by extrapolating the 
shallower models downward one or two points, reconverging, and then repeating the process 
until a sufficient depth in the atmosphere was obtained. Though the superadiabatic regions 
below the stellar surfaces are usually thin, a study of the development of convectively driven 
motions in flux tubes requires a substantial depth extent downward into the convection 
zone. For the sun, a depth extent of at least 5000 km is needed for the stability limit to be 
independent of mechanical conditions at the boundary locations (Rajaguru & Hasan 2000). 
For stars with log(gi) > 4.5, the ATLAS9 program fails to construct reliable models of sub- 
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surface regions upto depths that are required in our study. In those situations we attach 
adiabatic polytropic models that smoothly match with the models from ATLAS9; we note 
however that, in these cases, the realistic models from ATLAS9 cover the superadiabatic 
thin regions that drive the convective instability in the flux tubes. The further extension 
with a polytropic stratification is needed only to avoid breaking artificially the inertia of 
the convectively driven down fiow from the superadiabatic region above at the location 
where the shallower ATLAS9 model ends. The depth extents, and two important quantities 
characterizing the surface layers, viz., superadiabaticity 5 — V — Va^, where V and Vad 
are the actual and adiabatic temperature gradients respectively, and adiabatic index Fi, are 
tabulated in Table 1. The grid of models constructed have Te// in the range 4000 - 7000 K 
(in steps of 500 K) and log g in the range 2.5 - 5.0 (in steps of 0.5). 



2.2. Convective Stability Limits 

We have numerically solved the linear eigenvalue problem (equations [6] and [8] of Spruit 
& Zweibel (1979)), which determines the growth rates and frequencies of the unstable modes 
and their eigenf unctions, for stellar flux tubes using atmospheric models described in Sec- 
tion 2.1. The convective mode is the fastest growing and is identified from its cigcnfunction 
which has no node between the boundaries, corresponding to a monotonic downflow or up- 
flow throughout the extent of the tube. The critical values (3^ that determine stability against 
these convective modes are found. Assumption of temperature equality between the fiux tube 
and the ambient medium makes /? depth independent (Spruit & Zweibel 1979), and we eval- 
uate the critical field strengths from /3c using values of Pe at the height where the continuum 
Rosseland optical depth Tc—1. The above is a good approximation close to the photospheres, 
because the temperature differences between the fiux tube and ambient medium are minimal 
at these heights, as known from solar observations, and also the superadiabatic layers that 
drive the convective instability are very close to the photosphere (Table 1). The critical field 
strength s are the n given by Bc{g,Teff) = ^/87rpe/{l + Pc)= fcc(5',Te//)Beg(5',Te//). The factor 
fcc=l/ 1^(1 + f3c) represents the efficiency of CC; a value of fcc=l, i-e., /9c=0, corresponds to 
zero gas pressure, representing a fully evacuated fiux tube. The results are summarised in 
Figure 1. Shown are the critical field strengths Bc{tc = 1) (solid curve) required for con- 
vective stability and the maximum possible field strengths Beq{Tc — 1) (dashed curve), as a 
function of log g for various values of Teff- The variation of B^q with log g refiects the fact 
that the photospheric gas pressure, for hydrostatic equilibrium, varies as g. The dependence 
of Be on log g is modulated by the efficiency of CC, i.e., by icc{g,Teff)- Changing fee, i.e., 
varying efficiency of CC across the log g - Teff plane, would determine the trends of Be that 
differ from Beg in Fig. 1. To see this behaviour clearly, we have also shown the variation of 
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fee (dotted curves) with log g in Fig. 1. Several trends are seen in the figure. The most strik- 
ing feature is the increasing efficiency of CC for log g>3.5 as Teff increases; the maximum 
efficiency occurrs at log g=A and Te//=6500 K. For log g < 3.5 all stars show inefficient 
collapse, although 6 increases as log g decreases for a fixed Teff- Since gravity provides the 
main accelerating force on a displaced fluid element in a convectively unstable stratiflcation, 
a large reduction in gravity is expected to weaken the convective motions irrespective of 
variations in 5. This explains the inefficient CC at low values of log g. For log g> 3.5, CC 
shows varying amount of dependencies on Teff] the general feature is decreasing efficiency 
of CC with Teff- The largest differences between B^q and Be are found for the coolest star 
considered (Te^/=4000 K). For example, for Te//=4000 K and log 51=4.5, B^q and Be differ 
by about 1000 G. This inefficient behaviour of CC finds an explanation in Figure 2, where 
we plot the maximum values of the superadiabaticity 6 against Teff for different values of 
log g; 6 decreases monotonically with T^/f, for fixed log g, thus explaining the low efficiency 
of CC in cooler stars. As a general result, we find that CC is not efficient enough in stars 
cooler than Te//=5500 K, i.e., on the cooler side of the Sun in the main-sequence, owing 
to the absence of significantly superadiabatic sub-surface thermal structure, to produce field 
strengths close to Bgq. 

We emphasise that the stability limits found in this study are from a global analsyis 
and with effects of vertical gradients in ionization and adiabatic index taken into account. 
There are no known analytic criteria that capture the above effects. An analytic criterion 
derived by Roberts & Webb (1978) for finite extent vertical thin flux tubes is both necessary 
and sufficient but for linear (polytropic) temperature profiles and do not include effects of 
gradients in ionization and adiabatic index. Despite its limitations we use it to compare and 
understand quantitatively our results. This criterion yields critical values of P given by. 



where Hf, and Hq are the pressure scale-heights at the bottom and top ends of the ffux 
tube. We have compared values of f3c and i3c,rw in Table 2, adopting appropriate values 
for the quantities in Equation 1 from the stellar models used here (since V, and hence S, 
is assumed constant in the derivation of Equation 1, a nominal value of 2000 km is used 
as a representative extent of the flux tube). Though the exact values of Pe and /3e,RW do 
not match, owing to the reasons given above, they show similar trends against Teff. It 
is evident that the weaker held strengths (i.e., larger Pc,si) for stabihty in cooler stars is 
mainly a reflection of smaller 5. Additionally, since the constant temperature gradient of a 
polytrope is proportional to g, the scale-height ratio in Equation 1 yields small values (~ 1) 
at lower values of log g and hence large values of (3c- We attribute the differences between 
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the present results and those from Equation 1 to the neglect of gradients in V, ionization 
and adiabatic index in the derivation of Equation 1. It is noted that since magnetic pressure 
varies as and the gas pressure is proportional to difference between Be and Be^ would 
be larger in higher gravity stars for a constant value of (3c. 



3. Discussions and Conclusions 



We have examined the superadiabatic effect (or CC) (Parker 1978) in the stellar context, 
using realistic models of the outer convcctive layers of stars. Our results show that whereas it 
is possible to produce highly evacuated stable tubes in stars with Te// > 5500 K through the 
CC mechanism, it is not so in cooler stars; the decreasing amount of superadiabaticity in the 
upper convection zone layers of K and M spectral type stars make the CC inefficient yielding 
field strengths much less than Beg. Hence, if CC is the main physical process responsible 
for the formation of fibril-like evacuated flux tubes then it is expected that the differences 
between the observed B and Beg increase as Te// decreases. The existing observational 
results are not entirely consistent with such a trend as there are several cases of K and M 
dwarfs having B values close to or even exceeding Be^; but, all these stars with high B have 
high filling factors, /, typically 0.5 or higher, which correlate strongly with stellar angular 
rotation frequency Vt. Rotation-activity correlation dominated through high / values is a 
well observed fact (Saar 1990, 1996a) consistent with dynamo theory predictions. A careful 
look at the observational results compiled by Saar (1990) and Solanki (1992) reveal that all 
the cases of B < Beg also belong to K and M spectral types, but which have smaller Q and 
are less active; and all the G type stars show B values tightly around Beg. As an example, 
from the compilation by Saar (1990), we find that the stars GL 171. 2A (BY Dra) and HD 
201091, which are of similar spectral type (K5V, Te//Ri4400 K), have widely differing field 
strengths: 2.8 kG (Saar et al. 1987) and 1.2 kC (Marcy & Basri 1989) respectively. But 
these field strengths have good correlation with the Vl and / values (rotation periods of 1.85 
and 37.9 days and / values of 0.5 and 0.24, respectively). From our results shown in Figure 
1, we find for Te//~4400 K and a log gK,A.b (main-sequence) a field strength of ^ 1.3 kG 
that CC yields, in close agreement with the slow rotating, low / value case of HD 201091. 
Marcy & Basri (1989) themselves caution that the separation of B and / values is uncertain 
and that the flux /B is larger than expected for HD 201091. But, their speculation that 
B remains at the inferred value while / changes by 2 orders of magnitude, which is not 
inconsistent with the observed time variation in the chromospheric emission, flnds support 
from the present result. The best support for the present calculations comes from a multi- 
line infra-red Zeeman analysis of e Eridani, a K2 V star with Te//=5130 K and log 5f=4.7, 
by Valenti et al. (1995), who find /=0.088 and B=1.44 kG in close agreement with the value 
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Bc ~ 1.4 kG from Figure 1. It would thus seem that the CC indeed operates on all solar- 
like main sequence stars producing convectively stable tubes as dictated by the sub-surface 
superadiabatic structure; and, the higher than expected B values in K and M dwarfs may 
originate from effects induced by high values of Q. Recent detections of increased photometric 
variability in highly active K dwarfs, which exhibit saturation in their magnetic activity 
(O'Dell et al. 1995), indicate that stars with angular velocities Q > flgat show increased 
number of spots. Thus, the present result that Be -CBe^ for non-spot small scale magnetic 
fields in stars with Tg// <5000 K and its agreement with B observed on slow rotators provide 
a theoretical reason to believe the idea that flsat marks a change in contributions from / 
and B to /B (Saar 1996b): / saturates around 0.6 and Q(> ^sat) begins to contribute 
to B by increasing fspot/f to maintain the continued increase of flux /B well described by 
the power-law fits /Boc Q^-^ or /Boc (tcI^)^'^ to the observations (Saar 1990), where is 
the convcctive turn-over time and r^Q is the inverse Rossby number. We suggest that the 
calculated stability limits Be can be used to separate / from the observed fiuxcs /B for slow 
rotators {Q < ^sat), when there remain otherwise unknown uncertainties in the separation 
of / and B in observational analyses. The calculations here extend down only to Te/f—^OOO 
K due to limitations in generating reliable convection zone models for cooler stars, but, by 
extrapolation, it is likely that CC remains ineffective yielding Bc(tc — 1) <^ ^eq{Tc — 1) in 
cooler M dwarfs too. Two similar spectral type M dwarfs, observed by Johns-KruU & Valenti 
(1996), do show widely differing field strengths, in correlation with their Q, but the lower Q 
star seems to show B larger than the likely B^. Future refinements, both in observations and 
theory, are needed to understand fields in M dwarfs. 

We note that the thin tube approximation does not get constrained by increasing /; 
as long as the distribution of flux remains in the form of small flux tubes dominating the 
observed spectral profiles, the present results on B^ can be used in the interpretation of 
observations. Thus, if the observed profiles do arise only from non-spot fields even in fast 
rotating K and M dwarfs which show B^B^ then the present results call for new physical 
mechanisms to produce highly evacuated tubes, which are likely not due to thermodynamic 
reasons but due to fast-rotation induced phenomena. This may imply a failure of solar 
analogy as far as the formation and dynamics of surface magnetic fields are concerned. We 
conclude by pointing out that, in any case, the convective stability limits calculated here 
need to be satisfied by pressure confined flux tubes and thus may serve as useful lower limits. 

We are thankful to an anonymous referee for constructive comments. 
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Table 1: Characteristics of Model Atmosplieres 
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Table 2: Comparison of Critical Values of (3 

log g 

1.0 1.5 5 

Te// Pc [3c,RW PcRW Pc Pc,RW 

4000 1.81 7?76 2M 2^57 14.44 6.90 

4500 1.19 7.32 1.98 3.37 3.94 2.16 

5000 0.74 6.20 1.13 2.54 2.56 1.20 

5500 0.38 6.17 0.52 2.35 1.18 0.78 

6000 0.27 6.44 0.19 2.44 0.54 0.86 

6500 0.12 6.95 0.1 2.58 0.16 0.76 
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Fig. 1. — The critical magnetic field strengths Bc{tc = 1) for convective stability {solid 
lines) and thermal equipartition limits Beg(rc = 1) {dashed lines), against log g for various 
Tg//. Be correspond to those that can be achieved by convective collapse and Beg are the 
maximum possible field strengths for fiux tubes. Efficiencies of CC, fee, are shown as dotted 
curves and correspond to right-side ordinates of the panels. 
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Fig. 2. — The peak values of superadiabaticity 5=V-Va that characterize the sub-surface 
upper convective layers as a function of Teff ior three different values of log g. 



